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We have investigated how the dynamics of a quantum particle initially localized in the left boundary site 
under periodic driving can be manipulated via control of the right boundary site of a lattice array. Because of 
the adjustable coupling between the right boundary site and its nearest-neighbor, we can realize either coherent 
destruction of tunneling to coherent tunneling (CDT-CT) transition or coherent tunneling to coherent destruction 
of tunneling (CT-CDT) transition, by driving the right boundary site while keeping the left boundary site driven 
by a periodically oscillating field with a fixed driving parameter. In addition, we have also revealed that our 
proposed CDT-CT transition is robust against the second order coupling (SOC) between next-nearest-neighbor 
sites in three-site system, whereas localization can be significantly enhanced by SOC in four-site system. 

PACS numbers: 42.65.Wi, 42.25.Hz 


I. INTRODUCTION 

Coherent control of quantum dynamics via a periodically 
oscillating external field has been one of the subjects of long- 
lasting interest in diverse branches of physics and chem¬ 
istry [111 12] • One seminal result of the control is coherent de¬ 
struction of tunneling (CDT), a phenomenon originally dis¬ 
covered by Grossmann et al. in 1991 for a periodically driven 
double-well system ill], upon the occurrence of which the tun¬ 
neling can be brought to a standstill provided that the system 
parameters are carefully chosen. Due to its importance for un¬ 
derstanding many fundamental time-dependent processes and 
its potential application in quantum motor fllH] and quantum- 
information processing 0 ], CDT has received growing atten¬ 
tion from both theoretical and experimental studies. Theo¬ 
retically, it has been extended in various forms such as non¬ 
degenerate CDT 0L selective CDT isl], nonlinear CDT JH], 
many-body CDT ITol - lllll . instantaneous CDT lIT^ and mul¬ 
tiphoton CDT iflSn : Experimentally, it has been observed in 
many different physical systems like modulated optical cou¬ 
pler ifl^. d riven double-well potentials for single-article tun¬ 
neling iHa] . three-dimensional photonic lattices UM . a single 
electron spin in diamond |[l^ . Bose-Einstein condensates in 
shaken optical lattices ll2Qll2l|] . and chaotic microcavity j^ . 

In addition to conventional approach of modulating a lat¬ 
tice in a uniform fashion, modulating some certain lattice sites 
selectively also provides an attractive alternative for genera¬ 
tion of CDT, in which the rescaled tunneling amplitudes for 
different sites can be shut off effectively, thereby some in¬ 
triguing quantum manipulations are realizable such as dissi¬ 
pationless directed transport and beam splitter ll^ . CDT 
was traditionally thought to occur only at the isolated degen¬ 
eracy point of the quasi-energies and is related to dynamical 
localization i^ . However, recently a novel quantum phe¬ 
nomenon called dark CDT occurring over a wide range of 
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system parameters has been introduced in odd-A^-state sys¬ 
tems 11^ which is demonstrated to be caused by localized 
dark Eloquet state that has zero quasi-energy and negligible 
population at the intermediate state, rather than the superpo¬ 
sition of degenerate Eloquet states. Those advances on CDT 
studies offer benefits for coherent control of quantum dynam¬ 
ics. 

In this paper, we propose a scheme for the control of quan¬ 
tum tunneling in a periodically driven lattice array with con¬ 
trollable boundary through combination of characteristics of 
both normal CDT and dark CDT. We consider the coher¬ 
ent motion of a quantum particle in a lattice array, in which 
the harmonic oscillating external fields act on only the two 
boundary sites of chain. We have found an interesting re¬ 
sult that a single particle initially occupying the left bound¬ 
ary site experiences coherent destruction of tunneling to co¬ 
herent tunneling (CDT-CT) transition or coherent tunneling to 
coherent destruction of tunneling (CT-CDT) transition, when 
the driving amplitude of the external periodic field applied 
to the right boundary site is increased from zero. We have 
also revealed that the CDT-CT transition in three-site lattice 
is robust against the second order coupling (SOC) between 
next-nearest-neighbor sites, whereas the CT-CDT transition 
in four-site system is significantly affected by SOC. More¬ 
over, we present a good understanding of the results with help 
of the high-frequency approximation analysis as well as nu¬ 
merical computation of the corresponding Eloquet states and 
quasi-energies. 


II. MODEL SYSTEM 

We consider a single particle in an array of lattice sites with 
only two boundary sites driven by external periodic field, as 
shown in Eig. ([T]). In our model, the left boundary site is driven 
with fixed driving amplitude and driving frequency, while the 
right boundary site is driven with varied driving amplitude and 
fixed driving frequency. In the tight-binding approximation 
and assuming a coherent dynamics, the single-particle motion 
can be generally described by the tight-binding Hamiltonian 
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III. CDT CONTROL BY DRIVING THE RIGHT 
BOUNDARY SITE 


+ \N){N\ + 


N 


j=2 


N-l 


+ Y,^o(\j-l}{j+l\ + H.c.), 

j=2 


( 1 ) 


A. CDT control in three-site system 

We start our considerations from the three-site lattice array, 
the minimal one for odd-Wsite system. According to Eq. 
the evolution equations for the probability amplitudes read, 


where I 7 ) represents the Wannier state localized in the 
jth site, Qo is the coupling strength connecting nearest- 
neighboring sites, and vq is the second-order coupling (SOC) 
strength between next-nearest-neighboring sites. Instead of 
modulating a lattice in a uniform fashion, we modulate the 
on-site energies selectively. In this scheme, we assume a har¬ 
monic oscillating field applied in form of £’i (0 = Ai cos(ojt) 
for site 1, and EnO) = A 2 cos(ojt) for site N. Here Ai and 
A 2 are the driving amplitudes and oj is the driving frequency, 
respectively. 

We expand the quantum state of system ([T]) as |i/r(f)) = 
aj{l) \j), where aj(t) represents occupation probability 
amplitudes at the 7 th site, with the normalization condition 

2 = 1- From the Schrodinger equation10^(0) = 

j 

the evolution equation for the probability amplitudes 
aj(t) reads 


dai ^ ^ ^ 

i—— = El (t) a\ + i2o<22 + vo ^^3 
at 

(^7-1 + ^7+1) + ^0 (^7-2 + ^7+2) 

(7 = 2,3,...,A-1) 

da^ , . „ 

I—— = E^ (t) + yoaN -2 + ^oa^-i . 

dt 


( 2 ) 



FIG. 1: (Color online) Schematic of a lattice array with controllable 
boundary for model Q. Here the left boundary site is driven with 
fixed amplitude Ai and frequency cx), and the right boundary site 
driven with varied amplitude A 2 and fixed frequency o). All other 
sites are undriven. 


In this work, we mainly illustrate how the quantum dynam¬ 
ics of a single particle initially localized at the left boundary 
site can be controlled by driving the other boundary site, for 
both zero and non-zero SOC cases. 


= Ai cos (ojt) a\ + Qo <^2 + Vo <23 
= floai floas 

= A 2 cos (cot) as -h voai -h (3) 

We numerically solve the time-dependent Schrodinger 
equation © with Vo = 0 , i.e., neglecting next-nearest- 
neighbor tunneling in the chain. The initial state is (1,0,0)^, 
and the driving parameters of site 1 fixed as Ai = 22, co = 10 . 
The evolution of the probability distribution Pi = \ai\^ is pre¬ 
sented in figure (O for three typical driving conditions of site 
3. For A 2 / 6 C) = 0, Pi remains near unity, signaling suppression 
of tunneling. This is the quantum phenomenon well known as 
CDT. For A 2 / 6 C) = 2.0, Pi oscillates between 1 and ~ 0.4, 
showing partial suppression of tunneling. At A 2 /C 0 = 2.4, Pi 
oscillates between zero and one, demonstrating no suppres¬ 
sion of tunneling. The numerical results clearly indicate that 
the system undergoes a CDT-CT transition when A 2 /(o is in¬ 
creased from zero. Such a CDT-CT transition is more clearly 
demonstrated in Fig. ©(b), where the minimum value of Pi 
is used to measure the suppression of tunneling. When there 
is large suppression of tunneling, Min(Pi) is close to 1; when 
there is no suppression, Min(Pi) is zero. As clearly shown in 
Fig. ©(b), Min(Pi) slowly falls from its initial value to zero 
with A 2 IC 0 increasing from zero to 2.4. It is observed that the 
value of Min(Pi) drops to zero in narrow intervals around the 
zeros of /o(^ 2 /^), indicating that the particle is able to tunnel 
freely from site to site in these regimes. 

We explain the above numerical results through some 
analytic deduction based on high-frequency approxima¬ 
tion. In the high-frequency limit, we introduce the trans¬ 
formation bi = exp[-/JA i COS( 6 c) 0 <^^]< 2 i, Z ?2 = 

bs = exp[-/JA 2 cos( 6 c) 0 ^/^]^ 23 , where bj(t)(j = 1,2,3) 
are slowly varying functions. Using the expansion 
exp [ ±iA sin (cot) I co] = Y^kJk(^l <^) ( ±ikcot) in terms of 

Bessel functions and neglecting all orders except k = 0 in the 
high frequency region, we arrive at the effective equations of 
motion, 

= Do / o (^ i /<^)^2 

= ElQjQ(Ailco)bi + El{)J{)(A2lco)bs 

= kloJo(A 2 lco)b 2 . (4) 

Here we have dropped the SOC terms. It is easy to derive 
the analytical solutions of Eq. dU), 
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FIG. 2: (Color online) Three-site model 0]). (a) Time evolution of the 
probability at site 1, Pi = |aip, with different values of A 2 ; (b) The 
minimum value of Pi as a function of A 2 /du; (c) Numerical results 
of probability distribution Pi versus the driving amplitude A 2 and 
time, obtained from the original model 0 ; (d) Analytical results 
of probability distribution Pi versus the driving amplitude A 2 and 
time, given by the formula lH). The initial condition is {ai = 1,^2 = 
0,^3 = 01. The other parameters are chosen as Ai = 22, du = 10, 
^0 = l,yo = 0. 


b\ = -Cl + i — 1 ^=^= [C 2 sin (Kt) - C 3 cos (Kt)] 

Joi jj 2 , j 2 

Z ?2 = C 2 COS (Kt) + C 3 sin (Kt) 

b 3 =Ci- i—J^^ [C 2 sin i^Ki) - C 3 cos (Kt)], (5) 

//2 j 2 

where Cj(j = 1,2, 3) are constants to be determined by the 
initial states and normalization condition, /qi = Jo(M/oj), 

Jo 2 = Jo (M/oj) and K = -r Applying the ini¬ 

tial conditions /?i(0) = 1, ^ 72 ( 0 ) = 0, Z 73 ( 0 ) = 0 to Eq. ^ 
yields the undetermined constants Cj(j = 1,2, 3) in the forms 
Cl = C 2 = 0, C 3 = -i—=^L= and the occupation 

J 01+^02 V 4+4 

probability at site 1 as 


\bi\^ = 


42 +^4L^eos(7^0 


/2 + j2 

*^01 ^ '^02 


/2 + j2 

^01 ^ '^02 


( 6 ) 


From expression (| 6 ]), we immediately have two observa¬ 
tions: (i) when A 2 /tu 2.4,5.52, i.e., the zeros of / 02 , the 
minimal value of \bi^ is zero; (ii) when Ai/tu is tuned near to 
the zeros of /oi and A 2 /tu tuned far away from those values, 
|/?ip remains near unity. The analytical results of P\ based on 
formula ([ 6 l) are plotted in Fig. ©(d), which agrees well with 
the the numerical results obtained from the original model @ 
as shown in Fig. ©(c). 

The Floquet theory provides a powerful tool for under¬ 
standing the tunneling dynamics obtained in Fig. ©. Since 


Hamiltonian © is periodic in time, H(t + T) = H(t), where 
T = In I (jj is the period of the driving, the Floquet theorem 
allows us to write solutions of the Schrodinger equation © in 
the form aj(t) = aj(t) Qxp(-i 6 t). Here e is the quasi-energy, 
and aj(t) is the Floquet state. The Floquet states inherit the 
period of the Hamiltonian, and are eigenstates of the time evo¬ 
lution operator over one period of the driving 

U(T,0) = Tcxp[-i r H(t)dtl (1) 

Jo 

where 7“ is the time-ordering operator. Noticing that eigenval¬ 
ues of C(P, 0) are exp(-/ 6 P), the quasi-energies of this sys¬ 
tem can be numerically computed by direct diagonalization of 
U(T, 0 ). 

Our numerical results of the quasi-energies and the Floquet 
states for the modulated three-site system @ are plotted in 
Fig. @. There are three Floquet states with quasi-energies 
61 , 62 and 63 . We immediately notice from Fig. ©(a) that 
there exists a dark Floquet state with zero quasi-energy for all 
of the values of A 2 / 6 L). For this three-site system, there is no 
degeneracy in the quasi-energy levels. Therefore, the occur¬ 
rence of suppression of tunneling should be further explored 
through probe into the Floquet states. We display the time- 
averaged population probability (Fj) = (J^ dt\aj\^)IT for a 
given Floquet state (ai,a 2 ,a 3 )^ in Figs, ©(b)-(d). The Flo¬ 
quet state with (Pj) > 0.5 is generally regarded as a state lo¬ 
calized at the 7 th site. As seen in Figs, ©(c), the dark Floquet 
state has negligible population at site 2 while the population 
(Pi) > 0.5 holds for all values of A 2 I 0 J except those in the 
vicinity of zeros of Jq. The other two Floquet states are not 
localized at site 1 since their populations (Pi ) are lower than 
0.5. In a word, the CDT-CT transition shown in Fig. © comes 
from the dark Floquet state, whose population (Pi) undergoes 
a localization-delocalization transition. 



FIG. 3: (Color online) Quasi-energies and Floquet states of three- 
site system 0 . (a) Quasi-energies versus A 2 ICX). The time-averaged 
populations for the Floquet state in the quasi-energy level (b) e \, (c) 
62 and (d) 63 . The other parameters are Go = 1, vq = 0, u; = 10. 
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B. CDT control in four-site system 

We are now in the position to investigate the quantum con¬ 
trol of the four-site system. However, similar behaviors dis¬ 
cussed in next subsection will appear in six-site lattice system 
as well. For a four-site system, the coupled equations ^ read. 


dai ^ ^ ^ ^ 

i—— = Ai cos{ojt)ai i 2 o <22 + 
at 

da2 „ „ 

i -= ilQai i2o^23 -I- yoa4 

dt 

da^ „ „ 

i -= voai i2o( 22 + i2o<24 

dt 

da^ , ^ X ^ 

l—— = A 2 COS(6g)0^^4 + Vo< 22 + i2o<^3- 
dt 


( 8 ) 


We plot in Fig. ©(a) the minimum value of P\ versus 
by direct numerical simulations of the Schrodinger equation 
([8]) with Vo = 0. We start a particle at site 1 and fix the driv¬ 
ing parameters of site 1 as Ai = 11, co - 10. As shown in 
Fig. @(a), Min(Fi) takes extremely low values about zero 
except at a series of very sharp peaks. The quasi-energies 
of this system are shown in Fig. @(b), where we find that 
a pair of quasi-energies cross at the zeros of This is 

the origin of the extremely sharp peaks in localization seen in 
Fig. ©(a). Figs, ©(c)-(f) show that there is no localized Flo- 
quet states. These numerical results demonstrate the possibil¬ 
ity of inducing CT-CDT transition in four-site system through 
tuning the rescaled driving amplitude A^tio from zero to the 
points of quasi-energy crossing. 



FIG. 4: (Color online) Four-site model described by Eq. ©. (a) The 
minimum value of Pi as a function of A 2 /du. The extremely nar¬ 
row peaks appear at the zeros of J^iAijoS). (b) Quasi-energies versus 
Ail(o. At the zeros of JoiiA^joS), a pair of quasi-energies degener¬ 
ate. (c)-(f) The time-averaged populations (Py) of the Floquet states 
corresponding to 61 - 64 . The other parameters are Ai = 22, du = 10, 
Go = 1, Vo = 0. 


C. Effects of second-order coupling on CDT control 

In the above discussion, the infiuence of second order cou¬ 
pling (SOC), generally thought to be detrimental to CDT, is 
neglected. In this subsection. We have checked the robustness 
of our proposed scheme by direct numerical simulations of 
the Schrodinger equation © in the presence of SOC effects. 
As before, the left boundary site is initially occupied and its 
driving parameters is fixed as Ai = 22 , du = 10 . 

The influence of second order coupling on three-site sys¬ 
tem is plotted in Fig. ©. When SOC effects are taken into 
account, the numerical results in Figs, ©(a)-(b) show that the 
three-site system displays similar dynamical behaviors (CDT- 
CT transitions) as the case without considering SOC. The 
numerically computed quasi-energies and Floquet states are 
depicted in Figs, ©(c)-(f). Correspondingly, a localization- 
delocalization transition can be seen for the population dis¬ 
tribution (Pi) of the Floquet state corresponding to a quasi¬ 
energy 62 close to zero, when A^jco is increased from zero. 

We have also investigated the impact of SOC on the dy¬ 
namics of four-site system, as shown in Fig. ®. Fig. ©(a) 
shows the minimum value of Pi as a function of the driving 
amplitude A^jijJ. As A^IcjJ is increased from zero, Min(Pi) 
steadily decreases from a value of ~ 0.5 to zero before it peaks 
at A 2 / 6 G) 2.4. Contrary to our expectation, SOC facilitates 


rather than hindering the localization. This result is some¬ 
what counter-intuitive. From Fig. ©(b), it can be seen that a 
pair of the quasi-energies make a series of close approaches 
to each other as A^IcjJ increases. Detailed examination of the 
close approaches reveals that they are in fact avoided cross¬ 
ings. At the points of close approach, the tunneling is sup¬ 
pressed and the localization peaks. Compared with the zero 
SOC case shown in Fig. ©, the population distributions of 
Floquet states change dramatically for four-site system in the 
presence of SOC effects; see Figs, ©(c)-(f). The modification 
of the four-site-system’s dynamics due to SOC, such as en¬ 
hancement of localization, is the consequence of the Floquet 
states localized at site 1 instead of the quasi-energy degener¬ 
acy. 


D. CDT control in five- and six-site systems 

Our analysis above is given for three- and four-site systems, 
but similar behavior can be obtained for other finite number of 
sites when the SOC effects are not considered. 

The quantum dynamics of the driven A-site systems 
is investigated by direct integration of the time-dependent 
Schrodinger equation © with the particle initially localized 
at site 1. In what follows, we will present the numerical re- 




























5 



FIG. 5 : (Color online) Three-site model ([ 3 ]) with second order cou¬ 
pling. (a) The time evolution of Pi with different driving amplitude 
A2. (b) The minimum value of Pi as a function of (c) Quasi¬ 

energies versus AiIoj. (d)-(f) The time-averaged populations (Py) of 
the Floquet states corresponding to quasi-energies 61-63. The other 
parameters are Ai - 11 , oj- 10, Qq = 1 and vq = 0.2. 


suits for = 5 and 6 . Dynamics for five sites is depicted 
in Fig. o, which shows a similar CDT-CT transition as that 
of three-site system. Like the case of three-site system, this 
five-site system possesses a dark Floquet state with zero quasi¬ 
energy and negligible population at all of the even 7 th sites, as 
illustrated in Figs. (|7]) (c)-(d). The reason for CDT-CT tran¬ 
sition in the five-site system lies in the fact that population 
distribution (PQ for the dark Floquet state also experiences a 
localization-delocalization transition (see Fig. ©(d)). When 
N = 6, The CDT occurs only at the isolated points of param¬ 
eters as shown in Fig. ®(a), where a pair of quasi-energies 
become degeneracy (Fig. ®(b)). This is exactly the same as 
in four-site system. It is demonstrated that the CT-CDT tran¬ 
sition can be induced in six-site system by increasing A 2 /tu 
from zero. 


IV. CONCLUSION AND DISCUSSION 

In conclusion, we have studied how the dynamics of a sin¬ 
gle quantum particle initially localized in the left boundary 
site under periodic driving can be controlled by only driv¬ 
ing the right boundary site of a lattice array. In studying the 
number of lattice sites N = 3,4,5 and 6 , we have found that 
(i) a dark Floquet state with a zero quasi-energy exists in the 
three- and five-site systems. By raising the driving amplitude 


FIG. 6: (Color online) Four-site model (|8l) with second order cou¬ 
pling. (a) The minimum value of Pi as a function of A2ICX). (b) 
Quasi-energies versus A^joj. (c)-(f) The time-averaged populations 
(Py) of the Floquet states corresponding to quasi-energies 61-64. The 
other parameters are Ai = 22, u; = 10, Qq = 1 and vq = 0.2. 


of the external periodic field applied to the right boundary site 
from zero, we can induce a CDT-CT transition caused by a 
localization-delocalization transition of the dark Floquet state; 
(ii) no localized Floquet state occurs in the four- and six-site 
systems. However, we can realize a CT-CDT transition in the 
same way, thanks to the fact that a pair of quasi-energies de¬ 
generate at isolate points of parameters. We have also revealed 
that the CDT-CT transition in three-site lattice persists when 
the SOC effects are considered. Long-range interactions (non- 
nearest-neighbor couplings) in the lattice are important and 
non-negligible in some systems like biomolecules iH, poly- 
mer chains coupled w aveguid es m, and charge trans¬ 
port in a quantum dot array ||^[30]. Therefore, our proposed 
scheme provides a new route to possible application of CDT 
in such systems. Nevertheless, the CT-CDT transition in four- 
site system is significantly affected by SOC. In such a system, 
it is found, there exists a counter-intuitive phenomenon that 
SOC can enhance rather than hinder localization. 

Finally, we discuss the experimental possibility of observ¬ 
ing our theoretical predictions. The Hamiltonians ([T]) can 
be realized in different physical systems, for example using 
cold atoms or trapped ions in optical lattices elec¬ 

tron transport in quantum dot chains (lO, [H, [35|], and light 
propagation in an array of coupled optical waveguides with 
harmonic modulation of the refractive index of the selected 
waveguide along the propagation direction 1^ 1^. Recently, 
several z igzag-t ype waveguide systems were experimentally 
realized 11^13^ . which provides simpler realization of a one- 
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dimensional lattice with controllable first- and second-order 
(that is, next-nearest neighbors) couplings. Though experi¬ 
mental observation of CDT has been reported in many differ¬ 
ent physical systems, the proofs presented are not so rigorous 
as claimed, because CDT occurs only at isolated system pa¬ 
rameters. Since it is impossible to precisely determine a pa¬ 
rameter point through experiment, what is claimed is possibly 
not a genuine CDT which can be testified only after infinite 
time of evolution, something unrealizable in real experiment. 
Due to the limited time of evolution in experiment, the CDT as 
called may be a pseudo-CDT occurring around the surround¬ 
ing regime adjacent to an isolated parameter, instead of a real 
one. Hence, observation of CDT, such as the one in our con¬ 
sidered four-site and six-site systems, still present a generally 
unsolved challenge in current experimental setups 
Acknowledgments 


FIG. 7: (Color online) Five-site model, (a) Time evolution of the 
probability at site P\ - |aip, with different values of A 2 ; (b) The 
minimum value of Pi as a function of (c) Quasi-energies ver¬ 
sus (d) The time-averaged populations (Py) of the dark Floquet 
state corresponding to 63 . The other parameters are Ai = 22, u; = 10, 
Qo = 1, Vo = 0. 




FIG. 8 : (Color online) Six-site model, (a) The minimum value of Pi 
as a function of Ail(o\ (b) Quasi-energies versus Ailio. The other 
parameters are Ai -11, co- 10 , Do = 1 , vo = 0 . 
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